ODC  FILE  COP^ 


RADC-TR-80-1 00 

Interim  Report 
April  1980 


,  V-N 


/•A 

v  a 


4 


CO 


/vs 

W 


o 


Q 


ANALYTICAL  STUDIES  OF  PHASE 
ESTIMATION  TECHNIQUES 


Eikonix  Corporation 


Sponsored  by 

Defense  Advanced  Research  Projects  Agency  (DoD) 
ARPA  Order  No.  3223 


DTIC 

SELECT  Ef 

JUL  8  1980 


_  A 

APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


The  views  and  conclusions  contained  in  this  document  are  those  of 
the  authors  and  should  not  be  interpreted  as  necessarily  representing 
the  official  policies,  either  expressed  or  implied,  of  the  Defense 
Advanced  Research  Projects  Agency  or  the  U.  S,  Government, 


ROME  AIR  DEVELOPMENT  CENTER 

Air  Force  Systems  Command 

Griff iss  Air  Force  Base,  New  York  13441 


80  7  7  0  95 


This  report  has  been  reviewed  by  the  RADC  Public  Affairs  Office  (PA)  and 
is  releasable  to  the  National  Technical  Information  Service  (NTIS) .  At  NTIS 
it  will  be  releasable  to  the  general  public,  including  foreign  nations. 

RADC-TR-80-100  has  been  reviewed  and  is  approved  for  publication. 


APPROVED: 

DORIS  HAMILL,  1/Lt,  USAF 
Project  Engineer 


If  your  address  has  changed  or  if  you  wish  to  be  removed  from  the  RADC 
mailing  list,  or  if  the  addressee  is  no  longer  employed  by  your  organization* 
please  notify  RADC  (OCSE) ,  Griffiss  AFB  NY  13441.  This  will  assi.  us  in 
maintaining  a  current  mailing  list. 

Do  not  return  this  copy.  Retain  or  destroy. 


ANALYTICAL  STUDIES  OF  PHASE  ESTIMATION  TECHNIQUES 

Robert  A.  Gonsalves 
Robert  H.  Chidlav 


Contractor :  Eikonix  Corporation 
Contract  Number:  F30602-77-C-O176 
Effective  Date  of  Contract:  14  September  1977 
Contract  Expiration  Date:  14  December  1979 

Short  Title  of  Work:  Analytical  Studies  of 
phase  Techniques 
Program  Code  Number:  0E20 

Period  of  Work  Covered:  Jan  79  -  Aug  79 

Principal  Investigator:  Robert  A.  Gonsalves 
Phone:  617  273-0350 


Project  Engineer: 

Phone: 


1/Lt  Doris  Kanill 
315  330-3148 


Approved  for  public  release;  distribution  unlimited. 


This  research  vas  supported  by  the  Defense  Advanced 
Research  Projects  Agency,  of  the  Departnent  of 
Defense  and  was  monitored  by  1/Lt  ^is  Hanixl  (0CSE>. 
Criffiss  AF3  NY  13441  under  Contract  F30602-77-C-0176. 


j  Accession  i'er 

I  iris  Sfiil 

|  Di-J  *aa 

Cosax**  ■***•! 


S* _ 


V.^  -• 


i 


rj] 


UNCLASSIFIED 


•  rcnrocMiMS  o*SAi»it  at  ion  mmic  aho  aqowcss 

Eikonix  Corporation  s 

103  Terrace  Hall  Ave 
Burlington  HA  01803 


It.  CWI«OLUI|CO»riCt  MAAlf  MO  AODftttS  y 

Defense  Advanced  Research  Projects  Agenc/ // h  Apr®  1J80 
1400  Wilson  Blvd 
Arlington  VA  22209 


I*.  OtSTM»uYlOM  STATtwCMT  (*> 


Approved  for  public  release;  distribution  unlimited. 


17  C*lT*i»-jT>0«  ST  AT  Cm  CM  T  Cl  me  **t  Kmi  M.  II  mtUr—* 


I*.  VjmCMlHTART  MOTCS 

RADC  Project  Engineer:  Doris  Haoill.  1/Lt,  USAF  (OCSE; 


Phase  retrieval 
adaptive  optics 


•tlTRtCT  fC>  W»H  m  rtxrn  if  Mteeurp  mm  W»n»  If  wotwj 

Phase  retrieval  algorithns  are  presented  for  extended  objects  and  for 
nultiple  iaages  arising,  for  exacple,  froa  observations  at  several 
wavelengths  or  several  focal  planes.  The  first  experinental  deaonstra- 
tior.  of  phase  retrieval  is  reported.  ^ 


00  103  iaTio,on,«UrtB«wT( 


UNCLASSIFIED 

MCuoiTT  CvAUartCATrow  or  TmS Vaci  (*>«»  a«f»  Ewra^y. 


AQS 


UNCLASSIFIED 


KCu»itT  ciwuftCiTio*  o*  *— *  **6tnrk«  s«* 


1 


1  INTRODUCTION 

This  technical  report  summarises  the  work  performed  under 
contract  No.  F30502-77-C-0176  from  January  15,  1979  to  August  16, 
1979.  We  report  progress  in  the  following  areas: 

Draper  Work 

Algorithm  Development 

1.  SC  measurement  problem 

2.  Two  focal  planes  or  two  wavelengths  for  extended  objects 

3.  Smoothing  techniques 

Hughes  Worl^ 

1.  Noise 

2.  Blind  Test 

Experimental  Work 

One-D  Results 

2  DRAPER  WORK 

Personnel  from  Draper  gave  us  a  phase  aberration  that  is  a 
typical  output  from  a  simulation  of  the  ITEK  imaging  model.  Nc 
fit  this  aberration  with  a  14-term  Zemicke  polynomial  expansion 
and  constructed  a  point  spread  function  (PSF)  from  the  polynomial 
phase. 

When  we  performed  phase  retrieval  on  the  PSF  at  A  *  6  microns 
the  phase  was  retrieved  almost  exactly.  At  X  «  4  microns,  where 
the  effective  phase  aberration  is  twice  as  large  as  at  8  microns, 
the  algorithm  did  not  perform  as  well.  The  residual  phase  produced 
am  improved  PSF  but  it  was  far  from  diffraction  limited. 

At  am  August  16  meeting  of  Draper  and  EIKONIX  personnel,  we 
decided  to  perform  another  round  of  phase  retrieval  on  the  improved 
PSF. 

Details  of  this  work  appear  in  TM-9  and  a  report  of  the  successfu 
recovery  of  the  phase  in  the  above-mentioned  second  iteration  will 
appear  in  the  final  report. 


2 


3  ALGORITHM  DEVELOPMEN1 

3.1  DC  Measurement  Problem 

Quantization  of  the  observed  PSP  causes  an  error  in  the 
DC  (f  ■  0)  tern  of  the  optical  transfer  function  (OTP) .  In 
simulations  this  was  not  a  problem  because  the  data  was  generated 
in  floating-point  form.  However,  the  phenomenon  appears  in 
measurements  from  our  experiment  where  the  data  is  recorded 
only  to  3  significant  bits. 

The  current  algorithm  fits  an  observed  OTP  with  an 
hypothesized  OTP,  both  normalized  to  one  at  f  °  0.  Thus  an  error 
in  the  normalization  can  cause  an  incorrect,  hypothesized  OTP  to 
fit  the  observed  OTP  better  than  the  correct  OTP. 

The  remedy  is  simple.  Ke  find  the  constant  that  scales 
the  hypothesized  OTP  to  best  fit  the  observed  OTP  {in  a  mean  square 
error  sense) .  Calculation  of  the  constant  is  straightforward. 
Details  appeared  in  the  May  1979  monthly  report  and  will  be  updated 
for  the  final  report. 

3.2  Two  Pocal  Planes  or  Two  Wavelengths 

This  algorithm  is  particularly  useful  for  extended  objects. 
If  the  object's  spectrum  0(f)  is  known,  then  the  problem  is  nearly 
identical  to  that  of  the  point  object.  The  algorithm  uses  O(f)  and 
in  the  block  diagram  of  Figure  1.  P(f)  and  Pc(f)  are  the 
observed  and  hypothesized  OTF's. 

If  the  object  (scene)  is  unknown,  then  it  must  first  be 
estimated  based  on  additional  information.  Such  information  might 
cooe  fro®  an  additional  observation  of  the  scene  at  another  wave¬ 
length  or  et  isnuther  focal  position.  Call  the  two  observations 
x^tx)  and  with  Pcurier transforms  Sj(f)  and  22(f).  Then  an 

appropriate  nsci^ste  of  the  object  spectrum  is 


Pj^tfiz^f)  +  P2c*lf)z2Cf) 

JPlctf)|2  +  i*2ccf)|2 


where  P^=(f)  and  P2t.{f)  are  estimates  of  the  system  OTP's  based 
on  an  estimate  of  the  unknown  phase  aberration.  This  estimate 
is  used  in  a  slight  modification  to  the  algorithm  of  Figure  1  to 
find  the  best  phase  estimate. 

This  algorithm  was  tested,  successfully,  with  a  simulated, 
Gaussian  object  and  noiseless  observations  at  two  focal  positions. 
Testing  with  more  general  objects  will  be  reported  in  the  final 
report. 


Figure  1  Parameter  Search  Algorithm 


3.3  Smoothing  Techniques 

The  error  oetrio  of  Figure  1  has  many  local  minima  in  the 
N-dimensional  space  of  the  Zemicke  polynomial  coefficients.  In 
one-D  notation  the  metric  is 

£(C)  »  /|P(f)-  Pc(f)  |2  df, 

vhere  P{f)  is  the  observed  OTF  and  P^(f)  is  the  hypothesized  OTP 
with  parameter  vector  C.  In  earlier  work  under  this  contract  (TM-9) , 
we  developed  a  method  to  smooth  the  metric  by  averaging  over  C- 
intervals.  Thus, the  averaged  metric  is 

e7cT  =  /|p{f)|2df  -  2  Se{/P*Cf)P^THdf}  +  /|Pctf>|2df. 

where  the  overbar  denotes  the  C-average.  In  thrt  earlier  work, 
we  showed  plots  of  the  second  term  for  a  particular,  observed  OTP, 

P(f),  and  for  a  one-dimensional  variation  of  the  C-vector. 

In  this  reporting  period  we  studied  the  third  term  to  see  how 
rapidly  it  varies  with  the  several  C-paraneters.  TM-4  is  a  report 
of  this  study.  He  found  that  the  Zemicke  polynomials  for  coma, 

0  and  45  degree  astigmatism  and  X  and  Y  coma  were  relatively 
smooth  for  coefficients  from  about  .2  to  .6  (waves).  This  provided 
us  with  a  rationale  for  selecting  the  interval  (.25  to  .5)  for  the 
grid  search  of  the  basic  algorithm. 

The  result  above  is  a  by-product  of  the  search  for  a  good  smoothing 
technique.  He  have  neither  programmed  nor  tested  E(C),  itself,  to 
see  if  it  provides  a  better  metric  —  one  that  can  be  searched  more 
efficiently  because  it  has  fewer  minima.  This  task  has  not  been 
undertaken,  as  yet,  because  of  its  lower  priority  among  other  tasks 
to  be  accomplished. 


4  WORK  WITH  HUGHES  AIRCRAFT  COMPANY 


During  this  period  we  used  phase  retrieval  algorithms  developed 
under  this  contract  in  two  separate  tasks  for  Hughes  Aircraft  Company. 
The  first  was  a  study  of  the  algorithm’s  performance  with  additive 
noise  in  the  measured  PSF  and  the  second  was  a  blind  test.  At  the 
request  of  the  contract  monitor  and  with  Hughes'  permission,  ve 
briefed  the  government  on  this  work  on  August  16,  1979. 

In  the  first  task  we  estimated  the  residual  phase  error  vs. 
measurement  noise  in  the  PSF.  We  found  that  when  the  PSF  is 
corrupted  by  signal-dependent  noise  with  variance  kp(x),  it  causes 
the  residual  phase  (after  correction)  to  increase  gradually  up  to 
some  value  of  k  and  then  it  increases  rapidly.  Figure  2  shows  a 
typical  plot  of  the  residual  phase's  rms  value  vs.  input  noise  level 
when  the  original  phase  has  an  rms  value  of  0.3  waves.  This  plot 
was  obtained  by  Monte  Carlo  simulation.  For  k  higher  than  0.4  the 
residual  phase  increases  dramatically. 
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The  second  task  was  a  blind  test  where  Hughes  supplied  us  with 
a  PSF  from  a  system  with  unknown  (to  EIKONIX)  phase  aberration. 

Ne  estimated  the  phase  from  the  PSF.  The  results  are  summarized 
in  Figure  3.  Thu  real  phase  Cat  the  bottom)  was  given  to  EIKONIX 
by  Hughes  after  the  task  was  completed. 
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5  EXPERIMENTAL  WORK 

The  phase  retrieval  algorithm  has  been  extensively  tested  in 
computer  simulations,  in  which  a  point  spread  function  from  an 
aberrated  optical  system  is  generated  by  the  computer  and  the  PSF 
is  passed  to  the  phase  retrieval  program.  Although  this  kind  of 
work  is  a  necessary  first  step,  an  experimental  verification,  in 
which  an  optical  system  with  measured  aberration  is  used  to  create 
the  PSF,  is  needed  to  really  prove  the  concept.  An  experiment 
was  designed  to  provide  that  proof.  Ke  fabricated  a  phase  aberration, 
measured  it  interferoaetrically  (Figure  4) ,  put  it  in  an  optical 
system  to  generate  the  PSF  (Figure  5) ,  submitted  this  PSF  to  the 
phase  retrieval  algorithm,  and  compared  the  resulting  estimated 
phase  (Figure  6)  to  the  measured  phase.  The  agreement  is  excellent 
and  provides  the  first  actual  demonstration  of  phase  retrieval  in 
wavefront  sensing. 

A  detailed  presentation  of  the  experiment  will  be  given  in 
the  final  report. 


Figure  t  Interferoaetrically  measured  phase  aberration 
reconstructed  with  Zernicke  polynocial? 
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6  ONE-D  RESULTS 

Two  algorithms  were  developed  for  the  one-dimensional  phase 
retrieval  problem.  The  first  may  be  useful  for  small  phase 
aberrations  and  the  second  appears  to  be  a  satisfactory  solution' 
for  the  noiseless,  one-D  problem. 

Consider  a  small  phase  8(f)  for  which  the  coherent  system 
function  H (f )  has  the  approximation 

H{f )  =  eie(f)  ,  |f|  <  1. 

•  1  +  i8(f),  |f|  <  1. 

The  system  OTF  is  the  autocorrelation  of  E(f). 

More  generally,  we  define 

C(f)  =  1  +  ifl(f) 

where  0(f)  is  real  but  not  necessarily  small.  Its  autocorrelation 
P(f)  is  observed  and  the  problem  of  phase  retrieval  is  to  find  6(f) 

We  have  shown  that  this  problem  has  an  exact  solution  and  have 
developed  an  algorithm  to  find  "it.  The  solution  is 


K-l 

eK  ”  *N-K+1  “  iC  8  8  /6 

J»2  J  N+J-K 
X 

0N-K+1  °  BK-X+1  "  S  (6H+J-K  “  ' 

Ja2 

where  the  and  BR  are  samples  of  the  real  and  imaginary  parts  of 
the  measured  OTF  (less  the  diffraction-limited  OTF) ,  and  8^  is  a 
sample  of  8 (f ) . 

A derivation  of  the  result  and  several  examples  will  appear  in 
the  final  report. 
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The  second  algorithm  attacks  the  sore  general  problem  where 
the  sampled  OTP  is 


N-K+l 

E 

j=i 


b  *  a 

J  J+K-l 


We  can  show  that 


eN  -  FX  1  GK 


and 


0. 


K-K+l 


V 


where 


PR  “  angle  of 
Gk  "  cos-1  (Dj,/2) , 


and 


DK  ”  PN-K+1 


K-l 

E 


J»2 


u  “J+S-K- 


The  choice  of  signs  for  8_  yields  a  2-fold  ambiguity  for  each 
N  ** 

value  of  K  or  a  2  -fold  ambiguity  for  N  6^'s.  The  ambiguity  is 


Dk  has 


resolved  by  monitoring  the  magnitude  of  Dj..  If,  at  any  K, 
magnitude  greater  than  2,  an  incorrect  sign  has  been  made  prior  to 
the  current  value  of  K.  This  prunes  the  binary  tree  in  such  a  way 
that  the  number  of  branches  grows  approximately  linearly  with  K, 
an  entirely  manageable  search.  This  is  an  observation,  not  a 
theorem,  based  on  some  searches  of  size  N  «  10,  20  and  64. 

Details  and  examples  will  be  given  in  the  final  report. 
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